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LUNAR  IRRADIANCE  MEASUREMENTS 


BACKGROUND 


One  of  the  most  serious  problems  facing  the  Navy  today  is  that  of  trans¬ 
ferring  information  to  submarines  without  affecting  their  operational  capa¬ 
bilities.  Ideally,  such  a  communications  connectivity  should  be  available  at 
all  operating  speeds  and  depths,  and  should  not  recuire  any  external  receiver 
which  would  compromise  the  submarine's  strategic  or  tactical  function.  Until 
recently,  most  attempts  at  solving  this  problem  have  focused  on  systems  that 
use  the  lower  end  of  the  electromagnetic  spectrum  for  transmission  of  the 
desired  information  ( ie ,  elf,  vlf,  etc).  The  advent  of  high  power  ultra¬ 
violet-visible  lasers  has  chanqed  this  spectral  emphasis  to  include  optical 
fronuencies.  In  particular,  the  Navy  and  the  Defense  Advanced  Research 
Proiects  Agency  have  been  jointly  investigating  the  possibility  of  an  optical 

solution  for  strategic  communications.1 2 3 4  The  intent  has  been  to  assess  the 
practicality  and  suitability  of  optical  technology  to  solve  the  strategic 
communications  problem.  The  approach  involves  detailed  analyses  and  experi¬ 
mentation  in  the  areas  of  optical  engineering,  channel  characterizations,  and 


1  2 

critical  crmnonent  technology  development.  '  One  of  the  investigations  pur¬ 
sued  in  the  channel  characterization  portion  of  the  program  was  on  the  effect 
of  clouds  on  short-pulse,  optical  plane-wave  propagation  through  the  atmos- 
2,3,4 

nhere."'  '  This  report  describes  results  obtained  durinq  one  aspect  of  that 

investigation.  More  specif icallv,  a  partial  set  of  inferred  integrated 

optical  thickness  values  derived  from  lunar  irradiance  measurements  of  the 

laser- illuminated  cloud  layer  used  in  the  Downlink  Laser  Cloud  Propagation 

Experiment  performed  in  Kauai,  Hawaii,  during  Auqust  and  September  1979  is 

_  ,  2,3 
report  .<d . 

TECHNICAL  INTRODUCTION 


When  a  plane  wave  laser  pulse  is  incident  on  a  thick,  plane-paral lei 
cloud  layer,  the  resultant  pulse  is  severely  degraded  relative  to  the  incident 
illumination  profile  because  of  the  multiple  scattering  which  occurs  within 
the  cloud.  In  more  certain  terms,  an  initial  narrow  pulse  plane  wave  will 
exit  the  cloud  layer  decollimated  and  attenuated,  with  the  incident  pulse 


1.  Karp,  S,  A  Test  Plan  for  Determining  the  Feasibility  of  Optical  Satellite 
Communications  Through  Clouds  at  Visible  Frequencies,  NOSC  TN  279,  1  July 
107«. 

2.  Technical  Advisor  to  the  Strategic  Blue-Green  Optical  Communications 
Program  joint  Coordinating  Committee,  Naval  Blue-Green  Single  Pulse  Downlink 
Propagation  Model,  NOSC  TR  3R7 ,  Jan  1979. 

3.  Technical  Advisor  to  the  Strategic  Blue-Green  Optical  Communications 
Proqram  Joint  Coordinating  Committee,  Strategic  Blue-Green  Optical  Communica¬ 
tions  Program  Plan,  NOSC  Technical  Document  269,  16  July  1979. 

4.  Hostetter,  GR,  Test  Plan  May  1,  1979  Downlink  Laser  Cloud  Experiment, 
Contract  **N000 14-7B-C-07 1 6,  GTE  Sylvania,  Inc. 
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shape  qreatly  lengthened  and  disvcrted.  '  The  magnitude  of  these  effects 
depends  on  the  cloud  layer's  physical  thickness,  optical  thickness,  solar 
phase  function  and  rms  scatter  angle  (or,  conversely,  the  average  cosine  of 


2,7,8 


,8-18 


7,19,20 


seatterina) Several  analytical  and  Monte  Carlo  based  '  '  treat¬ 

ments  of  all  or  part  of  the  above  exist  in  the  literature.  Many  of  the  above 
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have  been  used  in  the  Navy  Single  Pulse  Downlink  Propagation  Model  to  simu¬ 
late  the  channel  effects  incurred  by  both  the  ground-based  and  space-base..', 
laser  communications  concepts  being  pursued  in  the  Strategic  Blue-Green 


5.  Stotts,  LB,  Closed  Form  Expression  for  Optical  Pulse  Broadening  in 
Multiple-Scattering  Media,  Applied  Optics,  17,  No.  4,504  (1978). 

6.  Dell-Imagine ,  RA,  A  Study  of  Multiple  Scattering  of  Optical  Radiation 
with  Application  to  Laser  Communications,  in  Advances  in  Communication 
Systems,  vol  II,  edited  by  VA  Balakrishnar  (Academic,  New  York,  1966).  p  1. 

7.  Danielson,  RE,  Moore,  DR,  and  van  de  Hulst,  HC,  The  Transfer  of  Visible 

Radiation  Through  Clouds,  J.  A^mosph .  Sci . ,  26,  1078  (1969). 

8.  Mooradian,  GC,  Geller,  M,  Stotts,  LB,  Stephens,  DH,  and  Krautwald,  RA, 
Blue-Green  Pulsed  Propagation  Through  Fog,  Applied  Optics,  18,  No.  4,  pp 
429-441  (1979). 

9.  Heqqestad,  HM ,  Optical  Communication  Through  Multiple  Scattering  Media, 
M.I.T.  Technical  Report  No.  472,  22  November  1978  (unpublished). 

10.  Bravo-Zhivotovskiy,  DM,  Dolin,  LS ,  Luchmin ,  AG,  and  Sarelyev,  VA,  Struc¬ 
ture  of  a  Narrow  Light  Beam  in • Sea  Water,  Atmospheric  and  Oceanic  Physics,  5, 
No.  2,  p  160-167  (1969)  (translated  by  PA  Kaehn ) . 

11.  Arnush,  D,  Underwater  Liqhc-Beam  Propagation  in  the  Small-Angle- 

Scattering  Approximation,  J.  Opt-  Soc .  Am.  62,  No.  9,  1109  (1972). 

12.  Fante,  RL,  Propagation  of  Electromagnetic  Waves  Through  Turbulent  Plasma 
Using  Transport  Theory,  IEEE  Trans  Antenna  Propagation,  AP-2 1 ,  pp  750-755 
(1973). 

13.  Ishimaru,  A,  and  Hong,  ST,  Multiple  Scattering  Effects  on  Coherent  Band¬ 
width  and  Pulse  Distortion  of  a  Wave  Propagating  in  a  Random  Distribution  of 
Particles,  Radio  Sci.  10,  637  (1975). 

14.  Honq,  ST,  and  Ishimaru,  A,  Two  Frequency  Mutual  Coherence  Function, 
Coherence  Bandwidth,  and  Coherence  Time  of  Millimeter  and  Optical  Waves  in 
Rain,  Fog,  and  Turbulence,  Radio  Sci.  11,  No.  6,501  (1976). 

15.  Stotts,  LB,  The  Radiance  Produced  by  Laser  Radiation  Traversing  a 
Particulate  Multiple-Scattering  Medium,  J.  Opt.  Soc.  Am.,  67,  No.  6,  p  815-19 
(1977). 

16.  Ishimaru,  A,  Theory  and  Application  of  Wave  Propagation  and  Scattering  in 
Random  Media,  Proceedings  of  IEEE,  65,  No.  7,  p  1030-61  (1977). 

17.  Fried,  DL,  Propagation  of  the  Mutual  Coherence  Function  for  an  Infinite 
Plane  Wave  Through  a  Turbid  Medium,  Optics  Letters,  1,  No.  3,  p  104-6  (1977). 

18.  Lutomirski,  RF,  Atmospheric  Degradation  of  Electro-Optical  System 
Performance,  Applied  Optics,  17,  No.  24,  p  3915-21  (1978). 

19.  Plass,  GN,  and  Kattawar,  GW,  Monte  Carlo  Calculations  of  Light  Scattering 
from  Clouds,  Applied  Optics,  7,  No.  3,  p  415-19  (1968). 

20.  Bucher,  EA,  Propagation  Models  for  Optical  Communication?  Through  Fog  and 
Clouds,  Proceedings  of  the  National  Electronic  Conference,  vol  29,  p  180-185 
(1974). 
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Optical  Communications  Program.  Unfortunately,  none  of  the  above  have  been 
experimentally  verified.  Hence,  a  technical  uncertainty  exists  in  any  perfor¬ 
mance  estimate  produced  in  the  system  engineering  portion  of  that  program. 

To  reduce  this  uncertainty,  the  program's  Joint  Coordinating  Committee 

initiated  an  experimental  effort"*'^  to  measure  the  transmission  and  multipath 

3 

time  spread  incurred  when  traversing  actual  atmospheric  clouds.  The  partici¬ 
pants  involved  in  the  experiment  were  GTE  Sylvania;  HSS,  Inc.;  and  the  Naval 
Ocean  Systems  Center  (NOSC) ,  with  NOSC  responsible  for  the  meteorological 
characterization  portion  of  the  experiment.  Because  of  the  importance  and 
magnitude  of  the  experiments,  NOSC  felt  that  an  alternative  characterization 
of  a  cloud  layer's  optical  thickness  was  desirable  and  that,  since  inhomo¬ 
geneities  in  the  cloud  layers  were  anticipated,  an  integrated  optical 
thickness  characterization  technique  should  be  employed. 

Gary  Lee,  et  al ,  of  McDonnell  Douglas  Corp.  suggested  that  a  method  using 
measured  solar  radiance  levels  could  be  used  to  infer  the  cloud's  integrated 
optical  thickness.  (See  reference  1,  appendix  B. )  This  procedure  is  general 
enough  to  apply  to  lunar  irradiance  as  well,  if  lunar  properties  such  as  moon 
phase,  and  earth-moon  distance  are  accounted  for.  The  author  was  tasked  to 
design  the  experiment,  perform  the  on-site  measurements,  reduce  the  data,  and 
document  the  findings.  This  report  provides  that  documentation. 


THEORY  OF  OPTICAL  THICKNESS  DETERMINATION 
USING  LUNAR  IRRADIANCE  MEASUREMENTS 

The  amount  of  radiant  energy  flux  from  the  moon,  as  measured  on  a  flat 
horizontal  surface,  depends  on  the  following  factors: 

21  22 

1.  Moon  phase;  a  full  moon  yields  the  most  light.  ' 

2.  The  angle  of  the  moon's  position  relative  to  the  zenith;  maximum 

21  22 

brightness  is  observed  when  the  moon  is  directly  overhead.  ' 


3. 


The  earth/moon  distance;  the  light  flux  is  the  most  intense  at  lunar 
22 

perigee . 


4.  The  matter  in  the  lunar  path  between  the  moon  and  the  detector,  such 
as  particulate  scattering  from  absorption  and  clouds.  Variation  of 
absorption  and  scattering  due  to  atmospheric  changes  is  to  be  con¬ 
stant  throughout  the  time  the  lunar  irradiance  measurements  are 
taken  . 


21.  Biberman ,  IM,  Dunkelman,  L,  Fickett,  ML,  and  Finke ,  RG,  Levels  of  Noctur¬ 
nal  Illumination,  Appendix  by  D.  R.  E.  Brown,  Institute  for  Defense  Analysis 
Research,  copy  109/200,  Contract  SO-SO,  Task  T-36. 

22.  Bond,  DS,  and  Henderson,  FP,  The  Conquest  of  Darkness,  AD  346297,  Defense 
Documentation  Center,  Alexandria,  VA,  July  1963. 
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Measurement  of  the  lunar  irradiance  with  a  clear  sky  and  measurement  of 
the  lunar  irradiance  with  clouds  in  the  light  path  can  be  used  to  determine 
the  scattering  and  absorption  due  to  clouds  if  the  moon's  phase,  the  moon's 
position,  and  earth/moon  distance  are  properly  accounted  for.  The  difference 
between  these  two  irradiance  measurements  is  the  essence  of  the  experiment. 

An  in  depth  discussion  of  the  theory  and  method  are  presented  in  reference  1. 

The  chief  advantage  of  this  method  of  determining  optical  thickness  is 
that  a  spatially  averaged  value  of  irradiance  is  obtained.  Thus,  if  the 
receiver  field  of  view  contained  some  patchy  clouds  and  some  clear  sky,  the 
value  of  irradiance  observed  would  be  approximately  equal  to  what  one  would 
see  if  the  visible  cloud  material  (of  sufficient  thickness)  had  been  equally 
distributed  over  the  entire  field  of  view.  This  will  not  work  for  clouds 
whose  optical  thickness  is  less  than  about  20,  because  the  method  relies  on 

7 

the  van  de  Hulst  equation  to  derive  r  from  a  transmission  value. 


MEASUREMENT  PROCEDURE 

Values  of  lunar  irradiance  were  obtained  utilizing  the  following  proce¬ 
dure  : 

The  irradiance  receiver  was  mounted  on  a  tripod,  exposed  to  the  lunar  sky 
with  a  zenith  orientation.  Values  of  lunar  irradiance  were  monitored  at 
various  times  throughout  the  night,  as  was  the  position  of  the  moon  at  the 
time  of  each  measurement.  The  values  obtained  were  used  in  conjunction  with 
the  theory  to  normalize  the  observed  irradiance  measurements  for  cloudy 
conditions  relative  to  those  taken  under  clear  sky  conditions. 

Simultaneous  with  the  lunar  irradiance  measurements  were  measurements  of 
direct  laser  pulse  propagation  through  clouds.  The  laser  was  mounted  on  an 
aircraft  and  pulsed  at  10  Hz.  The  aircraft  was  flown  above  the  clouds  and 
over  the  test  site.  The  position  of  the  aircraft  was  monitored  continuously 
throughout  the  aircraft  fly  bys. 

Figure  1  indicates  the  geometry  of  both  the  aircraft  laser  and  the  lunar 
irradiance  measurements.  The  laser  beam  was  spread  upon  exit  from  the  air¬ 
craft  and  further  dispersed  by  the  clouds  so  that  any  laser  light  indicating 
the  irradiance  receiver  was  negligible  and  not  detectable  above  the  moonliqht 
background . 

Measurements  were  taken  for  the  laser  experiment  during  the  last  week  of 
August  and  the  first  three  weeks  of  September  1979.  Lunar  irradiance  measure¬ 
ments  were  made  between  5  September  and  12  September  1979.  These  experiments 
were  conducted  simultaneously  on  the  Island  of  Kauai  on  Hawaii. 

EQUIPMENT  CONFIGURATION 

The  irradiance  receiver  consisted  of  an  RCA  4S17  photomultiplier  tube,  a 
high-voltage  power  supply,  and  a  preamplifier  buffer  mounted  in  a  waterproof 
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2  3 

housing.  In  addition,  a  45°  field  of  view  limiter,  a  diffuser,  and  a 
Wratten  filter  were  employed  on  the  receiver.  The  purpose  of  the  diffuser  was 
to  randomize  the  direction  of  the  incoming  light  rays  by  scattering  and  thus 
to  provide  a  spatially  averaged  effect.  Figure  2  shows  the  apparatus  along 
with  the  other  electronic  equipment. 


45 
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Figure  2.  Experimental  apparatus. 


The  45°  field  of  view  limitation  was  used  to  occlude  any  stray  light 
originating  from  city  street  lights  or  automotive  headlights.  Thus,  the 
measurement  was  not  one  of  irradiance  in  the  strict  sense  of  the  definition 
which  requires  a  180°  field  of  view.  Therefore,  a  factor  of  2  must  be  applied 
to  all  voltaqe  measurements  representing  cloud  attenuated  lunar  irradiance 
levels.  This  is  the  coefficient  which  appears  in  equation  5.  Details  of  the 
derivation  can  be  found  in  appendix  A. 

The  detector's  response  was  calibrated  over  the  entire  dynamic  range  of 
its  operation  usinq  a  standard  lamp  type  P-101  Electro  Optics  Associates,  and 
was  found  to  be  linear.  This  calibration  was  done  in  the  laboratory  with 
neutral  density  filters  of  varying  strengths  before  the  field  experiments  were 
begun . 

Data  acquisition  in  the  field  was  accomplished  by  reading  the  digital 
voltage  display  on  the  control  unit  and  simultaneously  storing  the  base  10  log 
using  a  strip  chart  recorder.  The  time  and  laser  position  (fig  1)  of  each 
measurement  were  noted  in  order  to  concurrently  correlate  the  results  with 
those  of  the  other  experiments.  The  moon's  angle  of  elevation  was  also 
monitored  to  the  nearest  degree. 


23.  Smith,  RC,  An  Underwater  Spectral  Irradiance  Collector,  J.  of  Marine 
Res.,  27,  341,  1969. 


DATA  REDUCTION 


The  method  of  obtaining  a  transmission  value  consisted  of  recording  the 
voltage,  (which  is  proportional  to  irradiance,  during  conditions  of  cloud 

cover),  with  the  axis  of  the  detector  pointed  at  the  zenith  (figure  1). 

A  standard  reference  voltage,  V  (1)  at  a  zenith  angle  of  ♦  was  recorded 

on  1  September  1979,  when  the  moon  was  observed  in  a  totally  clear  sky.  The 
presence  of  clouds  in  the  detector’s  field-of- view  during  a  calibration 
measurement  would  cause  scattered  light,  in  addition  to  the  direct  moonlight 
beam,  to  enter  the  detector  and  thus  would  lead  to  an  erroneously  high  value 
for  the  reference  irradiance. 


In  order  to  use  this  reference  value  to  calculate  transmissions  from 
irradiance  data  acquired  on  subsequent  nights,  it  is  necessary  to  correct  it 
for  chanqes  in  moon  phase,  moon  zenith  angle,  and  earth-moon  distance.  These 
corrections  are  summarized  in  equation  (1). 


V  (N)  =  V  (1  ) 
o  o 


P  (  N  )  Z  ( <|>  )  D(N) 
P(1)Z(*_)D(1) 

R 


( i ) 


In  equation  (1)  v^(N)  is  the  voltage  to  be  used  as  a  reference  on  night 

N,  where  N  is  the  day  in  the  month  of  September  in  this  case.  V  (N)  is 

proportional  to  the  irradiance  which  would  have  been  observed  from  the  moon  if 
the  sky  had  been  totally  free  of  clouds.  P(N)  is  a  factor  influencing  the 
irradiance  which  depends  upon  the  phase  of  the  moon  on  a  given  night,  N. 

P (N  )  =P (6  )=  1 . 0  for  a  full  moon,  which  in  this  case  occurred  on  9/6/79.  When  no 
moon  appears  in  the  sky,  P(N)=0.  On  the  other  nights,  the  lunar  irradiance 
will  be  a  fraction  of  that  of  the  full  moon,  measurements  of  which  are  de- 

2  1 

tailed  elsewhere  in  the  literature. 


Two  least,  square  polynomial  fits  to  these  data  are  given  by  equations  (2) 
and  (3). 


10 .  1787 

1.5813 

(N 

+  31  ) 

+  9.2956 

-2 

x  10 

(N  +  31 

)2 

-2.4818 

X 

10"3  (N 

3D3 

+  2.4894 

-5 

x  10 

(N  +  31  ) 

4_ 

(2  ) 

-2.6579 

X 

103  +  2 

.6784  x 

2  , 

10  (N  + 

31  )  -  1 

.0052  x 

1 

10 

(N  +  3  1  )  2 

+1.6733 

X 

10_1  (N 

+ 

3D3 

-  1.0416 

x  10"3 

(N  +  31  ) 

2 

(3  ) 

N  is  equal  to  the  day  in  the  month  of  September  when  a  given  experiment  was 
performed.  Equation  (2)  is  valid  when  the  moon  is  waxing,  in  this  case 
between  9/1/79  and  9/6/79;  whereas,  equation  (3)  applies  during  the  time 
period  between  9/7/79  and  9/12/79,  when  the  moon  was  waning. 


7,(<t>)  is  proportional  to  the  irradiance  measured  on  a  horizontal  surface 
at  a  qiven  zenith  anqle,  t .  Because,  to  a  very  qood  first  approximation,  the 
spectral  distribution  of  reflected  liqht  from  the  moon  is  equal  to  that 
emitted  from  the  sun,  solar  data  can  be  used  for  7,(i i> )  for  applications  where  a 

22 

ratio  7(  it  )  /?,(  ^  1  is  implemented.  7($)  is  qiven  by: 


loof/Un  =  2. BBSS  +  1.0P34  (i  -  $0  -  1.078b  (^  -  $  ) 2 


1.7092  X  10 


-  4>)3  +  1,2418  x  in' 


(4  ) 


and  is  most  accurate  for  zenith  anqles  between  0  and  47  decrees ,  with  4>  in 
units  of  radians.  The  data  for  which  this  expression  is  a  fit  were  also  taken 

2  1 

from  the  literature. 


D(N )  expresses  the  chance  in  illumination  resultina  from  the  variation  in 
the  eartb/moon  distance.  At  the  lunar  neriqee,  the  moonlioht  is  0.20  times 

no 

hriqhtor  than  at  the  apoqeo."" 

The  transm i ss i nns  were  obtained  by  normalizinq  the  corrected  voltaore 

readinqs  for  the  cloud  experiments  with  the  V  (N)  computed  from  equation  (2). 

o 


2Vi 
V  ( N  > 


o 


(S) 


where  T  -  transmission.  Values  of  t,  the  optical  thickness,  can  then  be  com¬ 
puted  usinn  an  equation  due  to  van  de  Hulst  and  later  improved  at  NOSC  to  be 
more  accurate  at  larqe  zenith  annles. 


A(  <t>  .  ) 

i 

0.  12fT 


11.10 


10  ) 


Here  <cos  f>>  ,  the  averaqe  cosine  of  the  sratterinq  anqle,  has  been  taken  to  be 
0.P7B.  <t>  is  the  zenith  anqle  of  the  moon,  in  radians.  (Bee  appendix  B.) 


'’14 

Aft.)  =  l.ftB  -  .^fll  <)>  .  +  2.7172  if  '  -  ft.ORftft  *  +  7.144ft  * 

i  j  i  l  i 


(7  ) 


S  ft 

-  3 . 4249  <t>  .  4-  P.  ft  1  SS  j>  . 

l  1 
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FINDINGS 


The  findinqs  are  in  accordance  with  the  precedinq  analysis.  Transmission 
values  are  qiven  in  tables  1  throuqh  S.  The  inferred  values  of  t  for  the 
reduced  measurement  data  are  listed  in  tables  6  throuqh  10. 

The  values  are  tabulated  by  GTE  Sylvania  run  number  and  position.  The 
run  numher  refers  to  the  value  assiqned  to  a  particular  overfliqht  of  the  air¬ 
craft  containinq  the  laser  and  is  used  as  standard  indicator  for  correlation 
between  the  various  experiments  conducted  at  Kauai.  Laser  position  is  that 
indicated  in  fiqure  1. 


No  reliable  data  were  obtained  durinq  conditions  of  rain  and  clear  sky, 
except  for  the  reference.  Some  run  numbers  are  therefore  omitted  when  there 
were  clouds  in  the  skv  surroundina  a  clear  path  between  the  moon  and  the 
detector.  The  data  for  these  particular  cases  are  not  valid,  for  reasons 
mentioned  in  the  previous  section  and  have  also  been  omitted. 


Values  of  T,  below  are  probably  not  very  accurate  because  those 

clouds  we  observed  durinq  the  indicated  time  period  were  too  thin  to  be 
adequately  described  by  the  van  de  Hul st  equation,  which  operates  in  the 
diffusion  reqime.  Al tbouob  this  is  the  case,  values  of  optical  thickness 
above  1  are  listed  for  the  purpose  of  rrrnparinci  these  results  to  numbers  for  T 
inferred  from  the  data  obtained  in  concurrent  experiments. 


The  wide  ranoe  of  values  obtained  for  the  optical  thicknesses  in  tables  G 
throuqh  10  reflects  the  wide  variety  of  clouds  observed  over  the  test  site. 
Gimilarlv,  the  variation  of  optical  thickness  within  a  sinqle  run  also  re- 
rip^ts  the  rapidly  channino  and  dynamic  behavior  of  the  clouds  we  encountered. 

Measurements  obtained  when  the  moon  was  full ,  or  near  full,  are  more 
reliable  than  those  acquired  on  niohts  with  a  lesser  deqroe  of  lunar 
i 1 1  urn inat ion . 


There  is  a  random  error  of  about  h9  on  each  value  of  the  optical  thick¬ 
ness.  The  back  around  radiation  from  the  ninbt  sky  and  stars  was  neqleeted,  as 
were  effects  orioinatinq  from  rbanqes  in  the  absorption  of  the  atmosphere  from 
dav  to  dav.  This  is  expected  to  be  a  fairly  oood  assumption  because  there  is 
not  much  atmosphere  absorption  in  the  spectral  reqion  of  530  nm .  Also  nealec- 
t ed  was  the  verv  small  chance  in  moon  phase  over  the  course  of  one  niqht's 
data  acquisition.  The  cosine  of  the  seat.terinq  anqle  was  assumed  to  take  the 
value  of  P.P75  (reference  7  and  appendix  B)  . 

RECOMMENDATIONS 

The  results  presented  in  tables  1  throuqh  10  are  for  only  one  method  of 
measurement  transmission  and  optical  thickness  of  clouds.  To  determine  if 
this  is  an  appropriate  method  of  measurino  inteqrated  optical  thickness  of 
clouds,  the  results  should  be  compared  with  other  measurements  and  postulated 
theory  for  cloud  optical  thickness.  Hence,  the  fol.lowinq  recommendations  are 
made  . 


1  1 


1.  Compare  optical  thickness  values  measured  by  this  method  to  those 
obtained  simultaneously  from  the  Knollenberg  particle  size  spectrometer  durinq 
the  blue-green  laser  cloud  propagation  experiments  performed  on  Kauai,  HI  from 

13  August  through  21  September  1979. 2,3,4 

2.  Correlate  the  multipath  time  spread  calculated  from  the  inferred 
optical  thicknesses  of  this  method  with  those  simultaneous  pulse  stretching 
results  to  be  contained  in  the  final  report  of  the  Downlink  Laser  Cloud 
Propagation  Experiment.  Several  models  for  multipath  time  spreading  (eg, 
references  2  and  5)  currently  exist  for  the  comparison. 

3.  Compare  pulse  shapes  predicted  by  Monte  Carlo  computer  simulation 
calculations  on  photons  transmitted  through  clouds  to  pulse  shapes  predicted 
by  the  present  values  of  optical  thickness  and  transmission. 
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7  500  yards  inbound  G.  7  500  yards  outbound 

5  000  yards  inbound  H.  10  000  yards  outbound 

J.  15  000  yards  outbound 
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5  000  yards  inbound  H.  10  000  yards  outbound 

J.  15  000  yards  outbound 


Table  4.  Cloud  Transmission  Values  from  Irradiance  Measurements 
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5  000  yards  inbound  H.  10  000  yards  outbound 

J.  15  000  yards  outbound 


Data  Acquisition  Period:  9/8/79,  23: 


Position 


32 


5  000  yards  inbound  H.  10  000  yards  outbound 

J.  15  000  yards  outbound 
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APPENDIX  A 


CORRECTION  FOR  USING  A  45°  ANGLE  FIELD  OF  VIEW  LIMITER  WITH  A  LAMBERTIAN 
RECEIVER. 

In  the  case  of  measurements  with  clouds,  the  irradiance  intensity 
incident  upon  the  horizontal  surface  of  the  collector  is  peaked  along  the 
vertical  direction  with  a  cos  9  dependence  for  the  intensity  of  light  rays 
whose  orientation  differs  from  the  zenith  by  an  angle  of  6. 


Figure  A-1 . 


The  total  amount  of  integrated  intensity,  1(0  )  received  at  the  collector 

L 


is  proportional  to: 

2w  e 


I  (0L)  a 


// 


0  0 


cos  9  sin  9  d9  d$ 


(A-1  ) 


The  solid  angle  is  accounted  for  by  the  integration  over  sin  9  d0  di)>, 

where  $  is  the  azimuthal  angle,  0  is  the  angle  of  the  field  of  view  limiter 

L 

when  measured  from  the  normal . 
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Note  that  it  is  not  necessary  to  normalize  the  irradiance  measurements 
because  a  ratio  between  I(ir/2)  and  I(ir/4)  would  cause  any  factors  common  to 
both  to  cancel  out.  However,  a  constant  of  proportionality  equal  to  1  A  is 
appropriate  in  order  for  the  integrated  irradiance  value  to  eoual  1  in  the 
case  when  no  field  of  view  limiter  is  present. 


2  IT 


I  (-=■)=- 

2  TT 


iff 


cos  0  sin  6  do  dij)  =  1  ( A-2  ) 

n  n 

The  cloud  is  assumed  to  subtend  the  field  of  view  in  this  analysis 

althouqh,  in  the  experimental  case,  the  irradiance  values  obtained  represent 

an  averaqe  of  the  sky  over  the  field  of  view. 

Next,  we  evaluate  the  inteqral  in  ( A  —  1 )  for  the  unobstructed  case,  with 

0  =  ir/2  and  then  for  the  experimental  case  when  0  =  ir/4. 

L  L 


I  (0  )  a  2tt 


i 
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cos  6  sin  ft  d0 


(  A-3  ) 
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cos  ft  sin  ft  d  ft  = 


2 

sin  9 


( A-4  ) 


There  fore : 


I  (j)  =  " 


(A-5  ) 


and 


I  (^-)  =  n.nnn  n 


(A-6) 


1  <i> 


1  '7> 


(A-7) 
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Thus  the  45°  field  of  view  limiter  is  seen  to  occlude  half  of  the  light 


transmitted  through  the  clouds  and,  as  a  result,  a  factor  of  2  is  obtained  for 
correcting  the  experimental  voltage  data  (which  are  proportional  to 
irradiance) . 
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APPENDIX  B 

LBS:sd 
Ser 81  1 ,24 
22  January  1474 

MEMORANDUM 

Erom:  Code  XI  14 

To:  SETCEile 

Via:  Code  81  I 

Subj:  Zenith  Angle  Dependence  of  Thick  Cloud  Transmittance 

Ref:  (a)  L.  B.  Stotts.  “Technical  Comments  on  the  Single  Pulse  Downlink  Propagation 

Model  Portion  at  the  Third  Interim  OSCAR  Report,"  NOSC  MEMO  Ser  XI  1  '560, 
dtd  2  January  074. 

( h )  R.  E.  Danielson.  D.  k.  Moore  and  IE  C.  van  de  Hulst,  "The  Transfer  of  Visible 
Radiation  Through  Clouds."  .1 .  Atmos.  Sci..  26.  No.  6.  pp.  1 078-X7  ( 1 464 ). 

( c >  E.  A.  Bucher.  "Computer  Simulation  of  Light  Pulse  Propagation  for  Communi¬ 
cations  Through  Thick  Clouds."  Applied  Optics.  12.  No.  10.  pp.  2341-2400 
(1473). 

(d)  E.  A.  Gray  hill.  An  Introduction  to  Linear  Statistical  Models.  McGraw-Hill.  New 
York.  1 46 1 .  Volume  1 . 

I .  The  purpose  of  this  memorandum  is  to  propose  a  new  functional  dependence  on  zenitn 
angle  for  thick  cloud  transmittance.  In  particular,  the  author  proposes  that  the  thick  cloud 
transmittance  in  the  Navy  Single  Downlink  Propagation  Model  he  given  by 

rc a  A (£c) .  ( 1 ) 

where 

6 

A<CS>  =  ^  ‘‘j  VSJ 

.1=0 


aQ  -  1 .64 

a,  =  -0.5513 
a2  =  2.7173 

a^  =  -6.4X66 
a4  =  7.1445 

a 5  =  -3.4244 
a(,  =  0.6155 

rather  than  by  tc  a  ( cos  “  (ref  (a)). 
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2.  The  zenith  angle  dependence  of  thick  cloud  transmittance  lias  been  modeled  by  Danielson 


Monte  Carlo  simulation. 

It  was  noted  in  reference  (a)  that  the  (cos<ps)^'"'  dependence  for  cloud  transmittance  now  used 
in  the  Navy  Downlink  Propagation  Model  breaks  down  for  zenith  angles  greater  than  84°. 

Using  a  triangular  factorization  (least-square  fit:  ref  (d))  of  the  X'X  matrix  composed  of  the 
tabular  results  given  in  reference  (b),  the  author  obtained  the  following  polynomial  expression 
for  the  zenith  angle  dependence  of  thick  cloud  transmittance: 


0 

A  <^s) 

=  ^  “j* 

j-o 

a0 

=  1  .69 

a| 

•T, 

»/•> 

d 

It 

U  *> 

=  2.7173 

a3 

=  -63)8(10 

a4 

=  7.1445 

a5 

=  -3.4249 

a6 

=  0.6155 

(3) 


Figure  It  I  is  a  plot  of  A(^  )  and  the  van  de  Hulst  Bucher  results  as  a  function  of  incident 
zenith  angle.  It  is  apparent  that  A(^rs )  is  a  good  approximation  of  the  van  de  Hulst  Bucher 
results  over  the  range  between  0  and  ‘■Hr.  It  is  recommended  that  this  expression  be  adopted 
in  the  Navy  Single  Pulse  Downlink  Propagation  Model. 


3.  Any  questions  on  the  above  should  be  directed  to  L.  B.  Stotts.  Code  SI  14.  X0858. 
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